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qq ■ Abstract 

The energy eigenvalues of the quantum particle constrained in a sur- 

■ face of the sphere of D dimensions embedded in a RP+ 1 space are 
i | obtained by using two different procedures: in the first, we derive the 

' Hamiltonian operator by squaring the expression of the momentum, 

■ written in cartesian components, which satisfies the Dirac brackets 
J> , between the canonical operators of this second class system. We use 
k> | the Weyl ordering prescription to construct the Hermitian operators. 

■ When D = 2 we verify that there is no constant parameter in the ex- 
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pression of the eigenvalues energy, a result that is in agreement with 
the fact that an extra term would change the level spacings in the 
hydrogen atom; in the second procedure it is adopted the non-abelian 
BFFT formalism to convert the second class constraints into first class 
ones. The non-abelian first class Hamiltonian operator is symmetrized 
by also using the Weyl ordering rule. We observe that their energy 
eigenvalues differ from a constant parameter when we compare with 
the second class system. Thus, a conversion of the D-dimensional 
sphere second class system for a first class one does not reproduce the 
same values. 
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1 Introduction 

When one obtains the representation of the physical operators given by the 
algebraic expressions of the Dirac commutators [|of a specific theory we must 
pay attention if these canonical operators are Hermitian ones. If it is not 
true, this fact indicates that there are the so-called operator ordering prob- 
lems. It is well known that many theories, for example, we can cite the non- 
linear sigma model fl|], the Skyrme model 0, D-dimensional sphere quantum 
mechanics f||, |j, which will be the subject of this paper, present ordering 
problems. The ordering problems, in many cases, appear as a consequence 
of the constraints dynamic acting in the systems. Therefore, the presence 
of the constraints can lead to a non-trivial relations for the Dirac commuta- 
tors which lead to an operator level problems for the representations of the 
physical operators of a peculiar theory. 

Some articles are devoted to study the problem of the quantization of the 
free particle constrained in a surface of D-dimensional sphere. We must men- 
tion the works of Abdalla and BanerjeeH who used a Lagrangian reduced 
method to obtain the quantum Hamiltonian of the multidimensional rotor; 
Kleinert and ShabanovQ] who employed the angular momentum algebra op- 
erators to calculate the energy spectrum; Foerster, Girotti and Kuhn|| who 
used the metrical formalism to described the quantum particle in a curved 
space. 

The purpose of this paper is to calculate with more rigour the energy 
1 These commutators define the quantum structure of a particular theory. 
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spectrum of a quantum particle lying in a surface of D-dimensional sphere. 
We use the second and first class Dirac methods of quantization^]. In both 
cases, the Weyl prescription is employed to construct the physical Hermitian 
operators. 

This paper is organized as follows. In Section 2, we derive the energy 
levels of the quantum Hamiltonian constructed from the Hermitian momen- 
tum operator, written in a rectangular coordinates, which satisfies the Dirac 
brackets. For D=2 we verify that this energy levels are in accordance with ex- 
perimental values[Q. This result can suggest the power of the Dirac method 
quantization of the constrained systems. In Section 3, we convert the second- 
class constraints into the first class ones. For this, we use the non-abelian 
BFFT formalism |7|], [Q . With the first class Hamiltonian we calculate the en- 
ergy eigenvalues of the quantum particle which differ by a constant parameter 
from the obtained by Dirac brackets. In Section 4, we give the conclusions. 
In the Appendices we perform a brief review about the non-abelian BFFT 
procedure, where we calculate the Lagrangian of this new first class system. 

2 Dirac brackets for the free particle on the 
D-dimensional sphere 

The quantization of the second class constrained systems is usually performed 
by using the method proposed by Dirac, Bergman and co- worker s||]. The 
constraints are classified as primary and secondary ones. Secondary con- 
straints are obtained from the condition that primary constraints are con- 
served in time. One must repeat the condition requiring time derivative of 
secondary constraints vanish until all independent constrains are obtained. If 
the whole second class constraints are established, the so called Dirac bracket 
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for the canonical variables A and B is given by 

{A, BY = {A, B} - {A, K}C~l{^ B}, (1) 

where Q and <pp are the second class constraints and the matrix elements 
C a p is defined by 

C af) = {<l>aM. (2) 

The quantum mechanics commutators are given by the replacement 
{,}*--*[,]■ 

Now, let us consider the dynamic of a particle in the D-sphere manifolds. 
The primary constraint is 

0i = XiXi - R 2 = 0, (3) 

where R is the radius of the sphere. With the expression of the classical 
Hamiltonian given by 

H c = ^PiPi , (4) 
we obtain the secondary constraint 

02 = XiPi = 0, (5) 

which expresses the fact that a motion on the surface of sphere has no radial 
component. From the expressions of primary and secondary constraints, (^) 
and d||), we obtain the algebra of the canonical operators which defines the 



quantum mechanics of a particle in a surface of D-dimensional spherepH,|12 
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[xi,Xj] = 0, (6) 

l 

[PjiPk] = -j^iPjXk-PkXj). (8) 

Here we would like to comment the operator ordering problem occurring in 
the right-hand side of Eq. (Q). This fact is solved under condition that this 
commutator satisfied the equation^ [pj,Pj] = . 

The solution of the momentum operator that satisfies the canonical rela- 
tions (0) and (H) is 



1 

Pk = - 



XkXi 



(9) 



The Laplacian operator in rectangular coordinates is [] 



Pk ■ Pk = d k d k + (OpOp + (D — l)Op) , (10) 

where the operator Op is defined as Op = Xidi and D is the dimension of the 
sphere embedded in a D + l dimensional cartesian space?]. The eigenstates are 
symmetric polynomials in the x^s (JTTJ defined as \polynomial) = jm^(x m + 



ix n ) 1 , where N(l) is a normalization factor, m and n take values m,n = 
1, 2, • • • D + 1, and I is an integer parameter^, / = 0, 1, 2, ... . The expression 
Pk ■ Pk is known as the laplacian V 2 on the D-sphere. 

2 We would like to acknowledge M. Plyushchay for the remarks about this question. 
3 The Laplacian expression in D + 1 dimension written in spherical coordinates is given 

by@ 

V 2 = £ ( E^l(nf =1 sin^)- 2 ( S in^)^ +2 - D ^( S in^^^) 

+ ags 1 sta*j)- a &). 

4 Notice that diXi = 8u = D + 1. 

5 The eigenvalues of the operator Op are defined by the following equation: 
Op\polynomial) = l\polynomial) . 
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We must mention the problem of ordering that appear in Eq. (0). Due 
to this fact, the canonical momentum pk is not a Hermitian operator as 
required by quantum mechanics. However, we can solve this problem by us- 
ing the Weyl ordering operator prescription 0. In the case of the canonical 
momentum, Eq. (H), this rule expresses that the new operator must be con- 
structed by counting all possible randomly order of the x's and d . Then, 
the symmetric momentum operator pk reads 



[pfc]sj/m (6<9/c j^x^Xidi j^x^diXi ^^XiX^di j^^idiXk 

diXfcXi ^ ®i"^i"^k) 



/ 



The quantum Hamiltonian is obtained performing the square of the Hermi- 
tian operator \p k ] S ym 



H 



[Pk]sym ■ [Pk] 



l -d k d k + ( OpOp + (D — l)Op 



D 2 -4 S 



2R 2 



(12) 



where the operator Op is defined in ref.flTDp. Thus, applying the Hamiltonian 
operator on the physical states \polynomial) we obtain the energy levels 



E, 



2R 2 



l{l + D-l) + 



D 2 -4 



(13) 



We would like to point out that the energy formula (|13]) gives no additional 
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constant energy for a particle on a circle^, D=2. This result is in agreement 
with the experimental energy level spacings of the hydrogen atom. Thus, 
for D=2, the Dirac quantization procedure together with the Weyl ordering 
prescription predict the correct energy levels for the quantum particle lying 
on the sphere. 

3 The non-abelian BFFT formalism for the 
free particle on the D-dimensional sphere 

In the appendix A we perform a brief review of the BFFT formalism and its 
non-abelian extension. In the appendix B we obtain the Lagrangian of this 
new theory. 

In the Section 2, we have seen the Hamiltonian for a particle of unit mass 
moving on the surface of a D-dimensional sphere of radius R isometrically 
embedded in R D+1 is given by 

H c = ^ViVi , (14) 
with the primary and the secondary constraints respectively written as 



(pi = XiXi — R 2 = 0, (15) 

02 = XiPi = 0. (16) 

The constraints T\ and T 2 are of the second class. The matrix elements of 
their Poisson brackets read 



6 For D=3 the extra term in Equation ( ^3|) is the same obtained in the three-sphere 
collective coordinates Skyrmions quantizationll3(]. 
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A Q/3 = {T a , T p } = -2e af3 XiXi , a, (3 = 1, 2 (17) 

where e Q/ g is the antisymmetric tensor normalized as ei2 = — e 12 = — 1. 

To implement the extended non-abelian BFFT formalism, we introduce 
auxiliary coordinates, one for each of the second class constraint. Let us 
generally denote them by n a , where a = 1,2, and consider that the Poisson 
algebra of these new coordinates is given by 

{r] a , rf} = u aP = 2e a/3 ; a = 1, 2. (18) 
From Eq. ( |A26| ), we have 

2X 11 X 22 = -2x i x i + Cl 2 T x . (19) 

After some attempts, we find that a convenient choice for these coefficients 
is 



Xu = R, 
X 22 = —R, 
X\ 2 = = A21, 
C\ 2 = 2, 

C\ 2 = 0. (20) 

Using (|53D , ([SEP, ( [Mil ), (0) and Q2IJ), the new set of constraints is found 
to be 



Ti = XiXi — R 2 + i??] 1 , 



(21) 
(22) 
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The first class constraint algebra is 

m,Ti} = o, 

{f 1 ,f 2 } = 2T 1 , 

{T 2 , T 2 } = 0. (23) 

Next, we derive the corresponding Hamiltonian in the extended phase 
space. The corrections for the canonical Hamiltonian are given by Eqs. ( |A2C| ) 
and ( A28 ). With the objective to simplify the expression of the first class 



Hamiltonian, we chose an algebra for the system defined by the parameters 
B b a in (|A27| ). We have verified that possible values are 



B\ = B\ = B\ = B\ = 0. (24) 
Using the inverse matrices 



U a /3 — l^-afii (25) 

x a e=(* (26) 

and the algebra defined by (|24D , it is possible to compute the involutive first 
class Hamiltonian 



H = \ P Ml - ^) - ^(1 - ^) + - (27) 

which satisfies the first class Poisson algebra 
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{T 1 ,H 2 } = 0, (B\ = Bl = 
{f 2 ,H 2 } = 0. {B\ = B\ 



0) 



0) 



(28) 
(29) 



Here we would like to remark that, contrary the results obtained by the 
abelian BFFT method applied to the non-linear Lagrangian theories [[l], |15|1 , 
the expression of the first class Hamiltonian (|27|) is a finite sum. As it 
was emphasized in the introduction, the possibility pointed out by Banerjee, 
Banerjee and Ghosh to obtain non-abelian first class theories leads to a more 
elegant and simplified Hamiltonian structure than usual abelian BFFT case. 

Here we intend to obtain the spectrum of the extended theory. We use 
the Dirac method of quantization for the first class constraints ||.The basic 
idea consists in imposing quantum mechanically the first class constraints 
as operator condition on the wave-functions as a way to obtain the physical 
subspace, i.e., 



T a \ip)p h y S = 0, a = 1,2. 



(30) 



The operators 7\ and T 2 are 



T\ = XiXi — R 2 + Rr) , 
T 2 = XiPi - Rrf + rfrf. 



(31) 
(32) 



Thus, the physical states that satisfy (|30D are 



\i)) P hys = tj S(xiPi - Rrf + r?V) 5(xiXi - R 2 + Rr] 1 ) \polynomial) , (33) 
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where V is the normalization factor and the ket polynomial is defined by 
\polynomial) = jmjfai -\-ix2) 1 ■ The corresponding quantum Hamiltonian of 
( JZTD ) will be indicated as 

Thus, in order to obtain the spectrum of the theory, we take the scalar prod- 
uct, p h ys {i})\H \ip) P hys , that is the mean value of the extended Hamiltonian. 
Then 



P h ys (ij\H\4>) phys = 

(polynomial] — J dn 1 dn 2 S(xiXi — R 2 + Rn 1 )5{xiPi — Rrf 2 + rfn 2 ) 

H5(xiPi — Rn 2 + n 1 n 2 )d~(x i x i — R 2 + Rn 1 ) \polynomial) . (35) 

Notice that due to the presence of the delta functions S(xiXi — R 2 + Rn 1 ) 
and S(xiPi — Rrf + rfn 2 ) in (|3~5D the scalar product can be simplified. Then, 
integrating over n 1 and n 2 we obtainQ 



phys(^\H\^) ph ys = 

(polynomiall^pXiXiPjPj - j^XiPiXjp ^polynomial) . (36) 

The final Hamiltonian operator inside the kets (|36"D must be Hermitian. Then, 

this Hamiltonian has to be symmetrized. Following the prescription of Weyl 

7 The regularization of delta function squared like (S(xiXi — R 2 + Rij 1 )) 2 and (5{xiPi — 
Ri] 2 + rfrj 2 )) 2 is performed by using the delta relation, (2tt) 2 S(Q) = lim^o J d 2 xe tk ' x — 
J d 2 x = V. Then, we use the parameter V as the normalization factor. 



J. Ananias Neto and W. Oliveira, 'Does the Weyl ordering prescription ...'12 



ordering (symmetrization procedure) we can write the symmetric Hamil- 
tonian as 



H S yrn — 2^> 2 \. XiX iPjPj]sym 2R 2 ^P^jPjlsym ' (^) 

where [xiXiPjPj] sym and [xiPiXjPj\ sym are defined as 



[xiXiPjPj] sym = ^{'iXiXiPjPj + AxiPjXiPj + AxiPjPjXi (38) 

+4p j x i x i p j + ApjXiPjXi + ApjPjXiXi), (39) 
1 

{xiPiXjPj] sym = — (2x i p i x j p j + 2x i pip j x j + 2x i p j p i x j (40) 

+2xiPjXjPi + 2xiXjPjPi + 2xiXjPiPj (41) 

+2p i x i x j p j + 2p i x i p j x j + 2p i p j x i x j (42) 

+2p i p j x j x i + 2p i x j p j x i + 2piXjXiPj). (43) 

Then, using the symmetric Hamiltonian operator H sym , eq. fl3"T|), the mean 
value (|35p is 
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phys\1p\H S y m \lp)phys 

1 1 

(polynomial] — [xiX iPjPj } sym - — [xiPiX jPj } sym {polynomial). 
1 

(polynomial\——r(AxiX i p j p j + AxiPjXiPj + Axip^Xi + ApjXiXiPj (44) 
4ort 

I Apj oc ^ j Ql/ ^ I 4^o ^ *^ ?■ *^ % / 

%P % X jP j I 2^2/ %P%P j X j | *2i3lj <iP jPi j I 23J iP j ^ jPi I 2*2'^ *^ jP jPi ^45 J 



48i? 2 



+2x i x j p i p j + 2p i x i x j p j + 2p i x i p j x j + 2p i p j x i x j + 2p i p j x j x i (46) 
+2piXjPjXi + 2piXjXiPj)\polynomial) . (47) 

The operator describes the momentum of free particle and its representa- 
tion on the collective coordinates space Xj is given by 



Substituting the expression 

as 



7r = —2 



g 

dXj 



(48) 



into fl44]), we obtain the energy levels, read 



E, 



l phys 



m 



sym | Y I phys 



2R 2 



l(l + D-l) + 



D(D + 1) 



(49) 



Comparing with the second-class energy expression, formula (0), we see that 
this expression differs by a constant parameter! 



8 Recently, a work of Hong, Kim and Parkjlj| suggest a generalized momenta definition 
that can conciliate the energy spectra obtained by Dirac and BFT formalisms, in the 
context of collective coordinates quantization of the Skyrme model. 
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4 Conclusions 

In this work we perform a study about the energy spectrum of a quantum 
particle constrained in a surface of the D-dimensional sphere. We use two 
procedures to obtain the D-sphere energy eigenvalues: in the first, we calcu- 
late the Hamiltonian by squaring the Hermitian momentum operator. This 
symmetrical momentum is derived, firstly, by an explicit solution of the Dirac 
brackets between the canonical operators, and then, employing the Weyl or- 
dering prescription, we obtain the Hermitian momentum operator. It is clear 
that the square of this operator (that is proportional to the Hamiltonian) 
also is Hermitian one. When we put D=2 (the physical system) the energy 
eigenvalues are according to experimental results. The second procedure is 
the conversion of the second class constraints of the D-sphere system into 
first class ones. We use an extension of the BFFT formalism. After this, 
we also employ the Weyl ordering rules to obtain the first class Hermitian 
Hamiltonian. Its eigenvalues differs by a constant parameter compared with 
the spectrum energy in the second class method. Thus, in the case of D- 
dimensional sphere quantum mechanics, Dirac Brackets quantization of the 
second class constraints can lead to the correct experimental results, contrary 
to the predict by the non-abelian BFFT formalism. 
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6 Appendix A: Brief review of the BFFT for- 
malism and its non-abelian extension 

The BFFT formalism of converting the second class constraints into the first 
class ones is a recent procedure?] used to quantize systems of second class 
constraints^. Then, the purpose of this appendix is to exhibit a review of 
the BFFT formalism. 

Let us consider a system described by a Hamiltonian H Q in a phase space 
{q % iP l ) with i — 1, . . . , N. It is supposed that there the system possesses only 
the second class constraints. Denoting them by T a , with a = 1, . . . , M < 2N, 
we arrive at the following algebra 

{T a , T b } = A ab , (Al) 

where det(A a b) ^ 0. 

As it was mentioned above, the general purpose of the BFFT formalism 
is to convert the second class constraints into the first class ones. This goal is 
achieved by introducing canonical variables, one for each of the second class 
constraint (the connection between the number of the second class constraints 
and the number of the new variables should be equal in order to preserve the 
same number of the physical degrees of freedom in the resulting extended 
theory). We denote these auxiliary variables by if and assume that they 
satisfy the following algebra 

K, V b } = " ab . (A2) 

9 Together with the Dirac first class quantization prescription. 
10 Compared with the well known Dirac brackets method. 
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Here u ab is a constant non-degenerate matrix ( det(u ab ) ^ ). The ob- 
tainment of u ab is embodied in the calculation of the resulting first class 
constraints which are denoted as T a . Of course, these constraints depend on 
the new variables r] a , that is 

f a = f a (q,p;rj), (A3) 
and are supposed to satisfy the boundary condition 

f a (q,p;0) = T a (q,p). (A4) 

In the framework of the BFFT formalism, the characteristic property of the 
new constraints is that they are assumed to be strongly involutive, i.e. 

{T a , f b } = 0. (A5) 
The solution of Eq. (|A5|) can be achieved by considering T a expanded as 

oo 

fa = Y.^ n \ (A6) 

n=0 

where is a term of order n in rj. The condition of compatibility with the 
boundary condition ( |A4|) requires 

= T a . (A7) 

Substituting the Eq. (|A^) into (|A5|) leads to a set of equations, one for each 
coefficient of rf 1 . We list some of them below 
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{Ta,T b } + {TW,T^} (rt) = (A8) 

{T a ,T b (1) } + {T«, Itf + {T^lf } (r?) + {Ti 2 ),T fe (1) } (r;) = (A9) 
{T a , if} + {Tj 1 ), T b (1) } ((?)p) + {Tf, T b } + {Ti 1 ), if > } (f?) 

+ {T?\T b {2) } iv) + {Ti 3 ),T 6 (1) } M = (A10) 



Here the notations {, }(q )P ) and {,}( v ), represent the parts of the Poisson 
bracket {, } corresponding to the variables (q,p) and (77), respectively. The 
equations above are used iteratively to obtain the corrections [n > 1). 
Equation ( |A5| ) gives Using this result together with the Eq. ([M]), one 
calculates T^ 2 \ and so on. Since is linear in 77 one can write it as 

T« =X ab (q,p) V \ (All) 
where X ab are some new quantities. Substituting this expression into ( TA8| ) 



and using (|A1|) and (IA4), we obtain 



A ab + X ac u cd X bd = 0. (A12) 

We notice that this equation does not define X ab in a unique way, because 
it also contains the still unknown elements uo ab . What is usually done is to 
choose u ab in such a way that the new variables are unconstrained. Con- 
sequently, the consistency of the method requires an introduction of other 
new variables in order to transform these constraints into the first class ones. 
This may lead to an endless process. It is important to emphasize that u ab 
can be fixed anyway. However, even if one fixes u ab , it is still not possible to 
obtain a unique solution for X ab . Let us check this point. A ab and uj ab are 
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antisymmetric quantities so expression ( |A12| ) includes M(M — l)/2 indepen- 
dent equations. On the other hand, since there is no additional symmetry 
involving X ab , they should represent a set of M 2 independent quantities. 

In the case when X ab does not depend on (q,p), it is easily seen that the 
expression T a + Tfp is already strongly involutive for any choice we make 
and we succeed in obtaining T a . If this is not so, the usual procedure is to 
introduce into Eq. (|A9|) in order to calculate and so on. At this point 
one faces a problem that has been the origin of some developments of the 
BFFT method, including the adoption of a non-abelian constraint algebra. 
This occurs because we do not know a priori what is the best choice we can 
make to go from one step to another. Sometimes it is possible to figure out a 
convenient choice for X ab in order to obtain a first class (abelian) constraint 
algebra at the first stage of the process M. It is opportune to mention that 
in ref . |18j , the use of a non-abelian algebra was in fact a way of avoiding to 
dealing with the higher orders of the iterative method. 

Another point of the usual BFFT formalism is that any dynamic function 
A(q,p) (for instance, the Hamiltonian) has also to be properly modified in 
order to be strongly involutive with the first class constraints T a . Denoting 
the modified quantity by A(q, p; 77), we then have 

{f a , A} = 0. (A13) 
In addition, A has to satisfy the boundary condition 

A(q,p;0)=A(q,p). (A14) 

The derivation of A is similar to what has been done in getting T a . Therefore, 
we consider an expansion of the form 
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(A15) 



n=0 



where is also a term of order n in 77's. Consequently, the compatibility 
with Eq. (|A14| ) requires that 



A™ = A. 



(A16) 



The combination of Eqs. (A6) 
equations 



([AT3D , ( |AT5D , and (1X161) gives the 



{r ai A} + {T a W A<% = (A17) 
{T Q ,A«} + {T a (1) , A} + {T«, + {Tf, = (A18) 

{T a , A< 2 )} + {T«, + } + {T«, A^} (v) 

+ {if, ^ (2) } ( „) + R (3) , = (A19) 



which correspond to the coefficients of the powers 0, 1, 2, etc. . . of the variable 
1]. It is just a matter of algebraic work to show that the general expression 
for A^ reads as 



A (n+1) = ^ ^ X bc Q{n) _ , A2Q n 

n + 1 

where uj a b and X ab are the inverses of uo ab and X a t, and 



n n—2 
m=0 m=0 



+ {T (n + l) ; (A21) 
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Finally, let us consider the case where the first class constraints form a 
non-abelian algebra, i.e. 



{T a , T b } = C c ab T c . (A22) 

The quantities C c ab are the structure constants of the non-abelian algebra. 
These constraints are considered to satisfy the same previous conditions given 
by (|E3D, QS3]), Q£§, and (|A7]). But now, instead of Eqs. M-flATD]), we 
obtain 



Cl b T c = {T^nj + iT^^}^ (A23) 
C c ab T^ = {T a ,T b ^} + {T^^T,} 

+ {r a (1) ,Tf } (?y) + {Ti 2 ),T b (1) } (r?) (A24) 

nc 71(2) _ fT". an(2)i 1 (rp(l) rr,(l)l 
{ ~'ab 1 c — l J a> J fc / + j J 6 J( g ,p) 

, r T (2) T (0h r T (l) T (3)| 

+{Ti 2) ,T b ( X + {Ti 3 ),T b (1) } w+ (A25) 



The use of these equations is the same as before, i.e., they shall work it- 
eratively. Equation (|A23|) gives T®. With this result and Eq. fl&24j) one 



calculates T^ 2 \ and so on. To calculate the first correction, we assume it is 
given by the same general expression ( |A11| ). Introducing it into ( |A23|) , we 
now get 

C c ab T c = A ab + X ac u cd X bd . (A26) 
Of course, the same difficulties concerning the solutions of Eq. ( A 1 2|) also 



apply here, with the additional problem of choosing the appropriate struc- 
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ture constants C c ah . To obtain the embedding Hamiltonian H(q,p,rj) one 
cannot use the simplified version discussed for the abelian case (embodied 
into Eq. ([A22[ ) ) because the algebra is not strong involutive anymore. Thus 
we start from the fact that the new Hamiltonian H and the new constraints 
T a satisfy the relation 

{f a ,H} = B b a f b , (A27) 

where the coefficients B b a are the structure constant of the non-abelian alge- 
bra. The involutive Hamiltonian is considered to satisfy the same conditions 
(|A14j )-( |A16|) . We then obtain that the general correction H^ n > is given by a 
relation similar to ( |A20|) , but now the quantities are given by 

n n—2 




m=0 m=0 

H^yM^-BjTW. (A28) 

7 Appendix B: the vacuum functional of the 
model and the Lagrangian that 
corresponds to the non-abelian first class 
Hamiltonian 

In this section, we intend to find the Lagrangian that leads to this new 
theory. A consistent way of doing this is by means of the path integral for- 
malism, where the Faddeev procedure |17| has to be used. Let us identify the 
new variables n a as a canonically conjugate pair (0, tt^) in the Hamiltonian 
formalism, 
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^-20, 

^ 2 -> 7I> , 



(Bl) 



satisfying ( |TH ) . Then, the general expression for the vacuum functional reads 



Z = N [djji\ exp{z / dt[iiPi + (fiir^ — H]} 



(B2) 



with the measure [do] given by 



[do] = [dxi][dpi][d(j)][d7r^]\det{,}\ 
5(xiXi -R 2 + 2R0)6(xiPi - Rtt^ + 2(jm^)J{5{K a ), 



(B3) 



a 



where A a are the gauge fixing conditions corresponding to the first class con- 
straints T a and the term \det{, }| represents the determinant of all constraints 
of the theory, including the gauge-fixing ones. The quantity N that appears 
in ( |B2|) is an usual normalization factor. Substituting the Hamiltonian ([27]) 
into (p2|), the vacuum functional reads 



Using the delta functions S(xiXi — R 2 + 2Rxj)) and S(xiPi — Rtt^ + 207^), and 
exponentiating the last one with Fourier variable £, we obtain 




(B4) 
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Z — N [ [d Xi ][d Pi ][d(j>][d^][d^\det{, }| S(x iXi - R 2 + 2Rxj>) J{5{K a ) 

J a 

exp{i J dt[xiPi - ip? (1 - ^) + 



+ 2* 1 - If) 



^(l-f)K)]}. (B5) 



Integrating over n^, we arrive at 



Z — N [dxi][d Pi ][d4>][dZ]\det{, }| <5(^ - R 2 + 2R(j))\{5{K a ) 

J a 

y^2$ exp{i J dt[~(l - ^)p 2 + (±i + x£)p. 

R 

1 00 R<j) 



2(l-f) 2 



(B6) 



Performing the integration over pj, we obtain 



Z — N [ [dx^d^Wdetl }\ 5(x iXi -R 2 + 2R<f>) Y[6(K a ) 



i-f\ 



1 



1 _^exp{i|^[ 2 - i _^ ) 2 (l-f)2 



+(l-^.)(x i x i + R^}}. (B7) 



Finally, the integration over £ leads to 



Z = N ( [dxi][d(f)]\det{, }\ 5(x iXi - R 2 + 2i?0)5(x i x i + J]^(Aa) 



! /" 1 1 (p(p /nn > 

i i? R \ R / 
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where the new S function that appear into the expression ( |B8| ) was obtained 
after integration over £. We notice that it does not represent any new restric- 
tion over the coordinates of the theory and leads to a consistency condition 



of constraint T\. From the vacuum functional (|B8| ), we identify the extended 
Lagrangian 

L= l - -1 % • (B9) 

2(l-f) 2(1-^)2 y 1 

Putting the extended variables, in the phase space, and equal to zero, 

we obtain the original Lagrangian. This result indicates the consistency of 

the theory. 
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